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Abstract. We analytically derive the expected number density distribution of 
Nambu-Goto cosmic string loops at any redshift soon after the time of string formation 
to today. Our approach is based on the Polchinski-Rocha model of loop formation 
from long strings which we adjust to fit numerical simulations and complement by 
a phenomenological modelling of gravitational backreaction. Cosmological evolution 
drives the loop distribution towards scaling on all length scales in both the radiation 
and matter era. Memory of any reasonable initial loop distribution in the radiation 
era is shown to be erased well before Big Bang Nucleosynthesis. In the matter era, the 
loop distribution reaches full scaling, up to some residual loops from the radiation era 
which may be present for extremely low string tension. Finally, the number density of 
loops below the gravitational cutoff is shown to be scale independent, proportional to 
a negative power of the string tension and insensitive to the details of the backreaction 
modelling. As an application, we show that the energy density parameter of loops 
today cannot exceed 10 -5 for currently allowed string tension values, while the loop 
number density cannot be less than 10 -6 per Mpc 3 . Our result should provide a more 
robust basis for studying the cosmological consequences of cosmic string loops. 



PACS numbers: 98.80. Cq, 98.70.Vc 
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1. Introduction 

It was shown by Kibble that cosmic strings are a natural outcome of cosmological phase 
transitions [TJ [2j [3], H] , and their possible presence in our Universe has been extensively 
studied ever since. Cosmic strings come in the form of line-like topological defects 
produced by the spontaneous breakdown of some symmetry [5j El [7j, but can also be 
stretched superstrings generated at the end of brane inflation [8| [9l [T0 | [TTl [12l [T3| [T^ | [15] . 
Once formed, cosmic strings evolve and very rapidly reach their so-called scaling regime 
in which the energy density of super-horizon sized strings behaves as radiation in the 
radiation era, or matter in the matter era [161 El EH EH] • Hence, for these long strings it 
holds that oc 1/d^, where d^ stands for the distance to the horizon (o?h oc t, the cosmic 
time). This scaling ensures that long strings never dominate the energy density of the 
Universe at late times and remain compatible with cosmological observations provided 
their energy density per unit length U is small enough. The current Cosmic Microwave 
Background (CMB) data limits the string contribution to at most 10% on the angular 
scales observed by the Wilkinson Microwave Anisotropies Probe (WMAP) satellite [20J, 
i.e. GU < 7 x 1CT 7 for Abelian vortices, G being the Newton constant [21] . However, 
being present all along the Universe's history, cosmic strings may be the dominant source 
of CMB non-gaussianities and are expected to dominate over inflationary perturbations 
at small angular scales [221 E31 EH ESJ [26] ET] EHl EHl ED]. 

The existence of a scaling regime for Nambu-Goto (NG) string networks requires 
the incessant formation of loops from long cosmic strings. Loops are produced at string 
autocommutation and intersection events, thereby evacuating the excess energy density 
from long strings. In the standard picture, the produced loops shrink and disappear by 
gravitational wave emission, which could provide additional observational evidence for 
cosmic strings [31j [32], [33], [3U [35]. Note, however, that Abelian Higgs string networks 
preferentially reach scaling by boson radiation rather than loop formation, at least 
during their numerically probed evolution [361 E3 138] . 

The expected cosmological distribution of loops has been subject to intense debate 
since the development of the original "one scale" model by Kibble [lj. In this model, 
at any time loops are formed with a size given by the typical distance separating long 
strings, i.e. a fraction of the distance to the horizon a\. NG numerical simulations in 
Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime have shown, however, that 
most loops are formed at much smaller size, suggesting the existence of tiny correlation 
lengths in cosmic string networks [TTJEHJEH]. These simulations are nevertheless limited 
in time and cannot exceed a redshift range typically of order 10 2 . Comparing this to the 
10 18 change in redshift expected at Big Bang Nucleosynthesis (BBN) [provided strings 
are formed at the Grand Unified Theory (GUT) energy scale], it is crucial for simulations 
to disentangle transient effects from the ones lasting over cosmological timescales. 

As shown in Refs. [IUJ [JT], many of the small loops observed in Nambu-Goto 
simulations indeed come from a transient effect associated with the relaxation of the 
initial string network towards its stable configuration. One must not expect these loops 
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to remain present after a long enough evolution. However, in addition to transient loops, 
Ref. [H] was the first to exhibit a sub-population of loops in scaling which constitute 
the cosmological attractor, a result recovered in Refs. [I2J H31 H3], albeit with some 
differences in the distribution. The energy density distribution of loops in scaling follows 
an evolution similar to long strings, i.e. it behaves as radiation in the radiation era, and 
as matter in the matter era. In this regime, the loop number density distribution has 
a power law shape: it is "scale-free" in the sense that loops of all size are present at 
any time. These results can be explained if the dominant mechanism of string evolution 
also applies to loops: they mainly loose energy by in turn emitting other loops. More 
quantitatively, denoting by a = i/d^ the size of a loop in units of d&, the loop number 
density n(a,t) and energy density p (a,t) in scaling read [H] 

dn _ S(a) dp _ jj ^ja) ,^ 

da a dl ' da a\ 

where the "scaling function" S (a) is well fitted, in the radiation and matter era, by the 
power laws S(a) = C Q a~ p with 




1 41 +008 
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Note that the energy density of loops evolves as does the one of long strings, p a oc l/d\. 

NG numerical simulations only account for string dynamics and intercommutation 
events, and in particular they know nothing about the impact of the loops' gravitational 
wave emission. In a realistic physical situation, a loop of size I emitting gravitational 
waves is expected to decay within a time roughly given by td = I /(TGU), where T is a 
numerical coefficient depending on the string structure, which for typical loops is of the 
order of T ~ 10 2 (13, [31]. As a result, for small loops with a < ad oc TGU, gravitational 
wave emission should overcome loop emission as the dominant energy loss mechanism. 
The (numerically found) scaling law of Eq. ([2]) can therefore no longer be used in the 
regime a < ad, and the loop distribution at those length scales must be understood 
further analytically, taking into account the emission of gravitational waves. 

Another phenomenon associated with gravitational radiation is the so-called 
backreaction effect which applies to both loops and long strings. The emission of 
gravitational waves renders the strings smoother and smoother on the smallest length 
scales, hence preventing self-intersection over a very small distance [161H7]. Therefore, 
long strings and loops cannot produce infinitely small loops: gravitational backreaction 
cuts off loop production below a certain scale. The typical length at which this happens 
was estimated to be of order £ c oc (GU) 1+2x t, where x > 0. When translated into 
the variable a and compared to the scale of gravitational decay ad oc TGU introduced 
above, this corresponds to a constant value of a c < ad [HI US [50] . 

The goal of this paper is to establish, both in the radiation and the matter 
dominated eras, a phenomenological description of the string loop number density 
evolution which includes the three above-mentioned effects: loop formation, decay 
by gravitational wave emission and the impact of gravitational backreaction. Various 
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analytical models have been proposed so far to describe the evolution of a cosmic string 
network [511 EH [531 EI] ■ Each of these models starts from a different set of assumptions 
and consequently, their results may differ. Our objective here is to describe the cosmic 
string loop distribution in the well-motivated situation where it coexists with a long 
string network in its scaling regime. Long strings indeed reach their cosmological 
attractor after a redshift range which does not exceed an order of magnitude [30] . For 
this reason, we use the Polchinski-Rocha (PR) approach of Refs. [551 EE] in which 
loops are formed from tangent vector correlations along long strings in scaling. Let us 
mention that a non-scaling loop production function has recently been introduced in 
Refs. [571 EH] to explain the existence of loops associated with the initial conditions. As 
discussed below (see Sec. I3.ip . our approach accommodates the existence of such loops 
without invoking a non-scaling loop production function. 

The PR model predictions are in good agreement with the correlators measured in 
Abelian Higgs simulations [38J and also with the NG loop distribution of Eq. (J2j) |59j, up 
to the overall normalisation. This agreement may appear surprising since simulations 
include many additional effects on top of loop formation from long strings, e.g. loop 
formation from loops, loop fragmentation and loop reconnection leading to daughter 
loops of a different size than the parent ones, all of them being non-negligible [18]. The 
physical interpretation is that these additional mechanisms preserve the functional form 
of the analytical PR loop production function (which is derived for the long strings) and 
amount to a renormalization of its amplitude. The PR production function can thus be 
considered as a robust physical feature against these effects, and for this reason, we take 
it as an ansatz for the total loop production function yet to be established. However, we 
can no longer fix its unknown coefficients analytically by using the long string scaling 
properties since we want to extend its use beyond the "long strings only" picture. We 
can, however, choose these coefficients to match the numerical results of Eq. (T5]), in the 
appropriate range of length scales. 

Below, we extend the method originally introduced by Rocha in Ref. [59] along 
the three following directions. Firstly, we introduce a change in behaviour of the loop 
production function below the gravitational backreaction length scale. The true shape 
of the loop production function in this regime being unknown, we consider various 
motivated possibilities and show that the final loop number density distribution remains 
unchanged. Secondly, we do not a priori neglect the transient solutions, which allows us 
to discuss the effect of the initial loop number density distribution at string formation 
time. This is the subject of Sec. [21 In Sec. El we discuss how the distribution evolves 
with redshift through the radiation and matter well as through the transition. 

We find that all loops rapidly reach the so-called "full-scaling regime", in which their 
number density assumes the universal form plotted in Fig. [2j We furthermore derive the 
loops' energy density parameter and the loop number density today before presenting 
our conclusions in Sec. HI 
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2. Cosmological attractor 

Borrowing notation from Ref. [59J, the distance to the horizon in terms of cosmic time 
t is expressed as 

= 3rb> ( 3 ) 

where the scale factor a(t) oc t v , with v — 1/2 in the radiation era and v = 2/3 in the 
matter era. Moreover, the cosmic time t can be expressed in terms of redshift z, of the 
Hubble parameter and of the density parameters today (H , Q T0 and fl mo , respectively) 
as 

^^i^iihf- (4) 

in the radiation era and as 

in the matter era. The above expressions constitute a very good approximation of the 
currently favoured ACDM (Lambda Cold Dark Matter) model [60J up to the smooth 
transition from the radiation to the matter era, and before the cosmological constant 
domination sets in (i.e. for redshifts z > 2). Note that, from Eqs. fll]) and fl5]), the 
instantaneous transition between the two eras occurs at 

~ " 1110 l=£(l + *q)-l, (6) 



16 ro 16 

where z eq stands for the redshift of equality between the radiation and matter energy 
densities. 

2.1. Evolution equation 

Denoting by n(£, t) the number density distribution of cosmic string loops of size £ at 
cosmic time t, we can write in the Eulerian description 

(T) 

where V(£,t) is the loop production function, i.e. the number density distribution of 
loops of size £ produced per unit of time at t. In the Lagrangian description, due to 
gravitational radiation, a loop of size £(t) shrinks at a constant rate such that 

dt - (8) 
which we take as the definition of 7d- According to the discussion of Sec. [fl and using 
Eqs. (j3J) and (jSJ), we have 

7d = ~ TGU. (9) 

Combining Eqs. ([7]) and (jHJ), one gets [59J 
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As suggested by the form of the scaling distribution in Eq. ((TJ), it is more convenient to 
express the loop sizes i in units of the horizon d^. In fact, since we shall be interested 
in the radiation matter transition, it is even more convenient to work in units of cosmic 
time t. For this purpose, let us define a new variable 7 and a new function J 7 , as 

f dn 

i(£,t) = -, Hi,t) = Tv (ii) 

respectively, in terms of which Eq. ffTUj) reads 

t ^gf 1 " (7 + 7d) = ahV^, t). (12) 

Evidently, the time variable t ranges from the initial time of string network formation 
(tim) to today (t ), while < 7 < 7 max , where 7 max corresponds to / max , the size of the 
largest (horizon-sized) loops, which will be discussed below. 

2.2. Loop production function 

The scaling law of Eq. ([[]) implies that, once the scaling regime is reached, t 4 J-"(7,t) 
should be a function of 7 only. From Eq. (fl2]l . we expect the same to happen for 
t 5 "P(7, t). Following Refs. [22 ESI ES], we moreover assume that this function is a power 
law, namely 

t 5 P( 7 ,t) = c 7 2 *- 3 , (13) 

where c and x are two parameters that will be fixed to fit Eq. ((2]). However, according to 
our discussion of Sec. [TJ the above expression is only valid for a range of 7 values where 
gravitational backreaction effects can be neglected. Hence, Eq. ( JT3l) holds for values of 
7 greater than 

7c = — — — ~ T(GU) 1+2 \ (14) 
1 — v 

where T is a number O(10) |50j . The parameter x i n this expression is expected to 
be the same as in Eq. ( 1131) . given that the gravitational backreaction length scale is 
precisely derived from the same two-point correlators used in the PR model to obtain 
Eq. ( !T3l) (see Ref. [50J). In our phenomenological approach, we consider 7 C as a free 
parameter below which the loop production function is given by a different power law, 

t 5 P(7<7c,t)=c c 7 2Xc " 3 . (15) 

The parameters c c and Xc are such that the loop production function is continuous 
across 7 = %, which leads to 

Cc = c7 c 2{x ~ Xc) . (16) 

Moreover, since gravitational backreaction smooths the strings on the smallest length 
scales, less loops should be produced in the small 7 limit, and one may naively expect the 
production function "P(7, t) to vanish for 7 = 0. In fact, by incorporating backreaction 
effects into the long string structure, it has been shown in Ref. [5D] that, for 7 ~ j c , 
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Eq. f)13p has a power law exponent renormalised to Xbr — 1 • This gives an extreme lower 
limit for \a 



as backreaction cannot smooth the string structure less than that. From Eq. ffl5|) we 
therefore see that the loop production function can be divergent in the 7 — > limit for 
1 < Xc < 3/2. However, being its first integral, the loop number density J r ( , y,t) will 
remain finite. Only in the limiting case Xc — 1 would we find a logarithmic 7—7-0 
divergence for J r ( , y,t) (see below), which nevertheless leads to a finite energy density 
(which is oc UEJ 7 ). This extreme situation corresponds to the production of an infinite 
number of infinitely small loops but with finite energy density. On the opposite, if 
gravitational backreaction was infinitely effective, no loops at all would be produced 
with 7 < 7c, corresponding to the limit Xc — > 00. Certainly the true physical situation 
lies somewhere between these two extremes [51J. In the following, we show that the 
final loop number density distribution is (almost) insensitive to the chosen value of Xc, 
and thus to the details of the backreaction smoothing process. Only the cutoff scale 
7 C enters into the expression for t 4 J-"(7, t) in its scaling regime. A sketch of the loop 
production function t 5 P(7,t) as a function of 7 is drawn in Fig. [TJ 

With Eqs. f|T3|) and (}T5|) we now have at our disposal a continuous expression for the 
loop production function for all 7 > 0. As mentioned earlier, however, there exists an 
upper bound 7 max on the accessible range of 7: any loop with 7 > 1/(1 — v) is of super- 
horizon size and hence indistinguishable from a long string. Hence, its contribution 
has already been included in the loop production function. In fact, in a cosmo logical 
string network, there does exist a small population of loops with size comparable to 
dh- We shall refer to them as "Kibble's loops", and they are not described by our 
current approach [HI H2]. Due to Hubble damping, their properties are in all aspects 
comparable to those of long strings, and they essentially decay by loop emission rather 
than gravitational radiation. As for long strings, we consider their net effect included 
in the loop production function. The typical size of these loops is given by the distance 
between two super-horizon sized strings in scaling, which is 



From the values of given in Ref. [UJ, one gets 700 ~ 0.32 in the radiation era and 
700 ~ 0.56 in the matter era. As we shall see later, our model predicts no loop larger than 
these values. However, in a realistic situation one should bare in mind that typically a 
few Kibble's loops are present. 

Finally, it has been shown in Ref. (61] that kinks propagating along the strings 
can also smooth the tangent vector correlator used in the PR model of loop formation. 
This amounts to setting x = 1/2 in Eq. f fl3|) for the 7 values at which this effect 
appears, i.e. it would require the addition of another piecewise domain in the definition 
of the loop production function. Although possible, we have chosen not to include kink 
smoothing in the following considerations since this effect is expected to be significant 



(17) 




(18) 
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in a domain 7 < 7k(£), where 7k (t) is a decreasing function of time [61] . As a result, the 
effect is transient and should no longer be observable when 7 k (t) hits the gravitational 
backreaction length scale j c . 

2. 3. Evolution from an initial distribution 

We now assume the loop number density distribution is known at some initial time £ ini , 
when it is given by 

AU£)^(£,t ini ). (19) 
at 

For instance, the initial loop distribution can be obtained by examining the statistical 
properties of the Higgs field phases just after the phase transition responsible for string 
formation [521 EH El] • However, since the loop production function of Eq. (|T5|) assumes 
a scaling form, A/ini(^) should be understood as the loop distribution expected at least 
just after the scaling of long strings takes place. In the following, we solve Eq. (TT2"j) . 
with the two- fold loop production function of Eqs. ( TT31) and f|T5|) . starting from a given 

A/ini at t = tini- 

Using new variables w and s defined from 

w = hxt, s = hit + ln(7 + 7d), (20) 
as well as Eq. (|T3l) for the loop production function, Eq. (fl2l) simplifies to 

A (a 3 ^) = caV 4 - (e— - 7d) 2x " 3 , (21) 

for 7 > 7 C . The same expression holds for 7 < 7 C up to the replacement c — > c c and 
X — > Xc- Assuming 7d constant and a oc t u , this equation is readily integrated in terms 
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of the Gauss hypergeometric function. For 7 > j c , the solution of Eq. (|2~T1) reads 

TV A c (7 + 7d) 2x ~ 3 „ „ 7d \ J( 7 t + 7dt) /00 ^ 
X(7, *) = Ti 2 F i 3 - 2x, //; /i + 1; — ■ + 3 , 22 

where X(x) is an arbitrary function (still to be determined) and 

M = 3i/-2 X -1. (23) 

For 7 < 7c, one obtains exactly the same expression up to the replacement of c and \ 
into c c and Xc [hence, fi c is defined as in Eq. ( 1231) by replacing X ~* Xc]- The unknown 
integration function in this regime is denoted by T c (x). In order to simplify the notations, 
we also define 

/(x)= 2 F 1 (3-2 X) /x;/x+l;a;), f c {x) = 2 F X (3 - 2 Xc , ^; Mc + 1; x) , (24) 
and the rescaled constants 

C = - , C c = - = JLtf<x-x») c (25 ) 

The loop number density distribution in its two respective domains then reads 

A ^ (7 + 7d) 2x " 3 J Id \ X[( 7 + 7d )t] 

X / 7>7c,t ) = C -4 / — — + 3 , 26 

t V7 + 7d/ a 

F^n, A r (7 + 7d) 2Xc ' 3 f / 7d \ X c [( 7 + 7d )t] 

^(7 < 7c, t) = C c - 3 ; c — ■ JH 3 . (27) 

r \7 + 7d/ a 

Continuity of ^(7, t) across 7 = 7 C at all times imposes that the unknown functions 
X[x) and I c (x) are not independent. From Eqs. (1261) and (l2Tj) one then obtains the 
relation 

• (7d + 7c) 4 



X c (x) =X(x) +K- 

x 



7c + 7d 



where K is a constant given by 



^ = C(7c + 7d) 2 *- 3 /(^H -C c (7c + 7d) 2Xc - 3 /c(^H. (29) 

Given J\fi n i(£), we can now uniquely determine the two unknown functions. At t — t ini , 
defining 

^ini = (7c + 7d)^ini , (30) 

one can use Eq. (T26]) to get 1{x) for all x > x^' 

X{x) = a 3 ni Mni(a: - 7dt ini ) - ( ' fX* "/(-tini)- (31) 

^ini V^ini/ V x ' 

This expression also provides I c (x) for all x > x; n i by using Eq. ( 1281) . Similarly, setting 
t = ti n i into Eq. ( 1271) uniquely determines X c (x) for all x < a;^, 

X C (X) = a^MAx - 7 dtini) " C c i± (f-) ^ /c(^tini) , (32) 

*ini VWni/ V X 7 



(2c 
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which now gives I(x < x ini ) from Eq. (1251) . After some algebra, plugging Eqs. (1251) . (1511) 
and (1321) into Eq. ( 1261) and (l2~Tj) . where the argument of the function X is now given by 
(7 + 7dK yields the final solution 



t A F{l>lc,t) 



UniJ \ a J 



7 + 7d I 1 - — 



+ C( 7 + 7d) 2x - 3 / 
-C(7 + 7d) 2x " 3 



7d 



7 + 7d 



t \ 2x+1 fa 



f 



7d ^in: 

7 + 7d t 



t 4 ^(7. <7<7c,t) 



a- \ 3 



t 4 -Af- ■ 

1111 ^ v mi 



7 + 7d 1 - 



+ c c ( 7 + 7d )^- 3 // 7d 

- C7( 7 + 7d) 2 *- 3 f- ) ( — ) / 



7c + 7d 
7 + 7d 



V7 + 7d 
t \ ' /a 

, ' 7 + 7d , 
a I 1 

7c + 7d 



7d h 



nu 

a / \7 + 7d * 

3 



a(t) 



t 4 J-(0< 7 <7r,t) 



t 4 ,hf- 



7 + 7d 1 - 



+ C c ( 7 + 7d) 2Xc - 3 /c 



7d 



7 + 7d 



C c (7 + 7d) 2 ^- 3 (-^ 



2Xc+l 



"1111 \ <■ 

— 7~ / J> 



Id t in 
7 + 7d * 



where we have defined 



lr{t) 



U 



(7c + 7d) — 



7d- 



(33) 



(34) 



Let us stress that, although the solution is defined on three different domains, it is 
continuous by construction. The "middle domain" 7 6 [7 r; 7c] owes its existence to two 
effects. Loops with a size I < j c t at time t can either come from "on-site production" 
with a 7 < 7c, or they may have started out earlier as larger loops, having 7 > j c , which 
then have shrunk into the domain 7 < 7 C by gravitational radiation. Of course, the 
same mechanism applies to all loops in the initial distribution, as it can be seen from 
the argument of Aim i n Eq. ( )33|) . weighted by a dilution factor a -3 . The third domain 
7 < 7 r corresponds to the virgin population of loops that started out with a 7 < 7 C and 
which knows nothing about shrunk loops produced at 7 > j c . As one expects, these 
"non-contaminated" loops cannot exist indefinitely and this domain of the solution only 
exists transiently. From Eq. (1341) . we find that j T (t T ) = for 



U - 1,. 
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Clearly, if 7 C <C 7^, this region exists only for a negligible duration after t mi . Even in 
the unlikely case of 7 C = 7^, one finds t T = 2t ini , i.e. at most one Hubble time. In the 
following, we neglect this 7 domain, i.e. we shall use only the first two solutions given 
in Eq. (133]) for the loop number density distribution. 



2.4- Asymptotic expansions 

We are now in a position to determine the unknown constants in Eq. (1331) by matching 
them to the scaling solutions obtained in NG numerical simulations |JT]. Since 
■A/i n i(£ > G?h|im) = 0, for any fixed 7 there exists a time at which all terms involving 
the initial loop distribution Af ini in Eq. (|33|) vanish. This relaxation from the initial 
loop distribution towards the scaling regime will be discussed in detail in Sec. El Taking 
the limit t ^> t ini and restricting our attention to the domain probed by numerical 
simulations, i.e. 7 3> 7d, one gets 

t^Hl > 7d, t > t ini ) ~ C 7 2x " 3 . (36) 
This expression matches with Eq. (pQ) provided that 

C = C (l-u) 3 ^, X = l-|- (37) 

As discussed in the introduction, the PR functional form of the loop production function 
therefore yields a scaling loop distribution which matches with the one observed in the 
numerical simulations. Fixing the value of the above coefficients to those given in Eq. (J2]) 
can be considered as a renormalisation procedure to include the additional sources of 
loop formation, such as fragmentation and reconnection, not considered in the genuine 
analytical PR model. 

In order to gain some intuition about the solution given in Eq. fl33|) . we also 
determine its asymptotic forms in the two following limits. On the one hand, at length 
scales smaller than 7^, but larger than 7 C we recover the result of Ref. |59j, 



Cu i 2x ' 2 

t 4 ^( 7c < 7 < 7d , t > t ini ) ~ 7T-7T i • ( 38 ) 

2 - 2x 7d 

Hence, in this domain the loop number density distribution is in scaling but with a 
power law exponent reduced by one. On the other hand, for loops smaller than the 
gravitational backreaction scale 7 C , we obtain for x c > 1 that 

t 4 J 7 (7<7c,t »tim) 



7c + 7d 



7d 

+ ^ r( ^ 2 +1)r fr, 2) . (») 

T(2xc -2 + fjic) 

which does no longer depend on 7. As a result, the loop number distribution under the 
gravitational backreaction scale also reaches a scaling regime and even becomes scale 
independent. Let us further mention the extreme case Xc — 1, which, as previously 
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discussed, yields a logarithmic divergence (in the radiation era, the matter epoch will 
be discussed in the next section), 



7d V 7 / 

The constant value of Eq. (139]) can be further expanded in the limit 7 C <C 7a- Using 
Eqs. (E5D and 022), one finally finds that 



which does no longer depend on Xc- Note that Eq. f HT]) smoothly connects onto Eq. ( |38|) 
at 7 = 7c, as it is to be expected. The loop number density distribution is sensitive 
only to the scale at which gravitational backreaction occurs (at 7 C ), but not how it 
occurs. This is fortunate since the precise form of the loop production function below 
the backreaction scale 7 C in unknown. For cosmological purposes, the only quantity 
that one needs to determine is thus 7 C . By comparing Eqs. ( l38i) and ( 14TI) . the physical 
interpretation is clear: the loop production function at 7 < 7 C is always much smaller 
than for 7 > 7 C (apart in the extreme case Xc — !)■ As a result, most loops of size 
7 < 7 C at any given time were originally produced with a 7 > 7 C and then shrunk down 
into the 7 < 7 C domain by gravitational radiation. When 7 C is not much smaller than 
7d, then \c encodes how smooth the change of behaviour from Eq. (|38|) to Eq. (HTj) is, 
but still does not affect the value of t 4 J-"(7 = 0,t) significantly (see Fig. [2]). 

In Fig. El we have plotted t 4 .F(7, £) in the radiation era and in the limit t 3> i ini , 
i.e. when neglecting all terms in Eq. fl33|) involving a dependence on the initial time t± ni . 
One recovers the three above mentioned domains with their power law behaviours of 
Eqs. (1361) . (1381) and (|4T|) . Note that the overall shape does no longer depend on cosmic 
time t but is universal: the loop distribution function is in its so-called "full scaling" 
regime. The black curve has been obtained with \ c = 2 whereas the dashed gray 
curves correspond, from bottom to top, to Xc = 4, % c = 1.5 and Xc = 1, respectively. As 
expected from Eq. ( HTj) . apart the logarithmic divergence when % c = 1, deviations due to 
the detailed shape of the loop production function below the gravitational backreaction 
scale remain negligible. 

Similar limits can be derived in the matter era by choosing C and x such that they 
match the relevant numbers of Eq. ([2]). The resulting loop distribution has the same 
overall shape as Fig. [2, with however different values for the power law indices. Notice 
also that the values of 7 C in the radiation and matter era do not need to be the same. 

In the next Section, we use the full solution given in Eq. (1331) to discuss how the 
loop number density distribution relaxes towards its attractors in both the radiation 
and the matter era. 

3. Relaxation towards scaling 

In this Section, we first discuss the relaxation of the loop number density from a given 
initial distribution A/iniW i n the radiation era. The motivation for these calculations 




(40) 



£ 4 .F(7<7c<7d,t>tim) 



2-2 X 7d 



(41) 
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Figure 2. The loop number density distribution t i ^F( i y,t) in the full scaling regime, 
i.e. when all terms involving i; n j in Eq. (|33l) can be neglected. The solution has been 
plotted for the radiation era (u — 1/2) with x an d C given by Eq. ([37]) and Xc = 2. 
The decay and gravitational backreaction scales have been fixed to the arbitrary values 
7d = 10 -4 and j c = 10~ 6 , respectively. The effects of changing Xc are illustrated by the 
dashed deviations which appear only for 7 < j c . From bottom to top, they correspond 
to the choices Xc = 4, Xc = 1-5 and the logarithmic divergent case Xc = 1- 

comes from the fact that just after the string formation time, the loop number density is 
expected to be different from its cosmological attractor. Therefore, we want to determine 
how long it takes to relax towards the universal form of Fig. |2j As we show below, for 
any reasonable initial distribution and string formation temperature, ^(7, t) is expected 
to be in full scaling well before BBN. The second problem concerns the relaxation of 
the radiation scaling loop distribution towards the matter era attractor. Put it in a 
different way, the question that one may ask is whether we can still observe traces of 
loops formed in the radiation era in the matter loop distribution. Our objective here is 
to express our results directly in units of interest in our observable Universe. For this 
reason, we express the cosmic time t in terms of the redshift z according to Eqs. (jlj) 
and 05]). We moreover use the currently favoured values of the cosmological parameters 
h = 0.72, n mo h 2 = 0.13 and Q V0 h 2 = 2.471 x 10" 5 [60], and all lengths and times are 
expressed in units of Mpc unless specified otherwise. Note that for these cosmological 
parameters we find from the definition of in Eq. ([6]) that z* ~ 3059. 
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3.1. Radiation era 

As a well-motivated example, let us assume that the initial loop number density 
distribution is given by the Vachaspati-Vilenkin (VV) random walk model of Ref. |65j, 

4^niW = a(^) 5/2 . (42) 

Here, the constant C\ depends on the initial correlation length £ of the string network, 
an order of magnitude estimate being given by [65] 

From Kibble's argument, causality imposes that £/dh(£ini) < 1 an d G\ > 1. A maximum 
value for C y can be estimated by imposing that the energy density of the string network 
at formation saturates the radiation energy density, i.e. 

H 2 (t ini )^—-^. (44) 
Neglecting constant coefficients, one then gets 

-^-) ~GU, (45) 



tim 



such that 



1 < Q < (GU)- 3/A . (46) 

The loop distribution of Eq. (j4"2j) makes sense only for £ > £, otherwise no loop can be 
produced. Also £ < d^ ini , with dh|ini = 2tmi in the radiation era (since the long strings 
are already included in the loop production function, see the discussion above), and 
A/ini(^ > 2ti n i) = 0. The maximum value of t A J\f- in i then occurs for £ = £ and can be 
estimated by 

tUfUt) * cf /3 < {cuy\ (47) 

By examining Eq. ( |33|) , let us first remark that since M- m \ is maximal at a fixed physical 
length £ — £, any effect associated with the initial loop distribution should be peaked at 
a non-constant 7f(i) = £/t and then travel towards smaller and smaller 7 values. This 
is what is observed in numerical simulations during relaxation of the initial conditions 
towards the loop scaling regime |JT] • The above derivation does not consider the damped 
evolution epoch that is supposed to take place just after the phase transition responsible 
for string formation [3]. During this period, the background radiation is expected to 
smooth strings and loops such that some of the small loops should disappear. The above 
estimate should therefore be viewed as a upper bound for A/Li- 

Using the VV initial distribution of Eq. (I42p . it is immediate to check that, at fixed 
7, the time dependence of all terms involving AAni in Eq. ( |33l) cancels. In other words, 
the initial loop density distribution remains, during some amount of time, on an equal 
footing with the freshly produced loops. However, since fif^i > 2t ini ) = 0, this regime 
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can only be transient since the vanishing initial distribution at super-horizon loop sizes 
will "sweep leftwards" towards smaller 7 as time goes on. For any fixed value of 7, all 
terms involving A/jni in Eq. ( |33|) vanish after a time th given by 

m = i±3i. (48) 

tini 7 + 7d 

The relaxation time th/t- m i is therefore a function of the scale of interest 7. For 7 ^> 7d 
(or 7d = 0), it varies as 2/7 such that it takes more time for the smaller loops to reach the 
scaling attractor, another property observed in the numerical simulations of Ref. |^Tj . 
For the smallest values 7 <C 7d, the relaxation time saturates to th/^ini — 2/7d> which 
is completely determined by the gravitational decay length scale. In terms of redshift, 
Eq. (SHI) reads 



1 + ^(7) _ /7 + 7d (4Q) 



1 + z ini V 2 + 7d 
which is limited from below by 

— > VlU? • (50) 

^ini 

Since 7d can be small, it may take quite a long redshift range for the initial loop 
distribution to vanish and the observability of such transient effects also depends on 

^ini • 

If one assumes that the correlation length £ is approximately £ ~ 1/yU, then 
Eq. OS]) implies that m p; t ini can vary from l/y/GU to 1/(GU), with m Pi denoting the 
Planck mass. Taking the latter value and using Eq. gives the lower bound 



1 fm 



>^H vm ^ lvm} - (51) 



Combined with Eqs. © and ( 1501) . one gets 



* > l^f-w ~ l ° 32GU - (52) 

For string tensions of current cosmo logical interest, GU > 10~ 10 such that > 10 22 : 
the decay of the initial loops is over well before BBN (z nuc — 10 10 ). In fact, for all string 
tensions such that GU > 10~ 28 (y/U > 100 TeV), it is safe to claim that any initial loop 
remnants have been radiated away at the time of equality between radiation and matter 
(z eq w 10 4 ). Using the lowest value t ini = £, the same statement would apply for strings 
as light as GU ~ lO" 37 , i.e. VU ~ 1 GeV. 

We have just shown that all terms involving Af ini in Eq. fl33l) can be reasonably 
neglected at all redshifts z < Zh- All the other terms involving t; n i exhibit a power law 
decrease, up to a correction factor by the hypergeometric function which rapidly equals 
unity. In the radiation era, the biggest terms in Eq. (|33|) are those appearing in the 
expression for 7 > 7 C and evolve as 

^ \ 2 X +1 . _ . , / \ 1-4y 
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Figure 3. Relaxation towards scaling of the loop number density in the radiation 
era from a Vachaspati-Vilenkin configuration at Zi n i = 10 19 . The step, whose position 
sweeps leftwards over the entire range of 7 with decreasing redshift, corresponds to 
the distance to the horizon at the string formation time above which A/ini = 0. It 
disappears after a finite time given by Eq. (|48| such that it is no longer visible at 
z = 6x 10 16 and the loop number density distribution is on the cosmological attractor. 



Using the value 1 — 4% ~ 0.2 from Eq. (j2J), they contribute less than 10% as soon as 
z/zmi < 10~ 5 , i.e. for 

z < z % = 10 26 VGU, (54) 

using Eq. (loTj) . Again, for all relevant values of GU, these terms cancel well before BBN. 

We have plotted in Fig. [3] the numerical solution stemming from Eq. (1331) for an 
initial loop distribution given by Eq. ( l4~2l . For illustrative purposes, we have chosen 
GU ~ 10~ 14 such that the string formation redshift is at least z ini ~ 10 19 while letting 
the gravitational effect at the same value as before. The relaxation towards scaling is 
completed at z ~ 10 16 , as one expects from Eq. f}54")) . 



3.2. Matter era 

It follows from the previous discussion that the cosmic string loop number density 
distribution expected at the end of the radiation era is given by its scaling solution 
described by the power laws of Eqs. ( |36|) . ( |38l) and ( 14TI) . Taking this shape as the initial 
distribution AC at z = z* [see Eq. (jSJ)] for the subsequent evolution in the matter era, we 
can now discuss how the loop number density relaxes towards its new scaling solution 
in the matter era. 

For v = 2/3, it follows from Eq. fTSBl that jj, — 1 — 2% and the hypergeometric 
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function in Eqs. f !24p simplifies to the polynomial expression 

/( *> = (T^T (> - 2^) ' < 55 > 

and similarly for f c (x), provided Xc 7^ 1- In the matter era, the extreme case Xc = 1 
would produce a logarithmic divergence in time for 7 < 7 C suggesting that it is no 
longer physically acceptable. Therefore, this extreme case will not be considered in the 
following. Concerning the values of 7d and 7 C , one does not expect the former to change 
significantly from the radiation to the matter era [as one can see from Eq. (jUj)]. However, 
given that 7 C is closely related to the tangent correlator [see Eq. (fl4|) ]. it assumes a lower 
value in the matter dominated epoch. 

As for the radiation era, it suffices to discuss the relaxation of the loop number 
density distribution in the regime 7 > 7 C since the smaller structures relax even faster. 
We denote by x R , C R and x M > C u the respective values of \ an d C in the radiation and 
matter era. The time dependence associated with the terms involving TV* (the matter era 
equivalent of A/lni for the radiation epoch) in Eq. (1331) shows that, unlike when starting 
from the VV distribution, they are now damped during matter domination and become 
sub-dominant compared to the scaling solution as soon as z < z^y), with 

1 -U y /Vi (r \ V(3/2-3 Xr ) 

1 + ^7J \ ^ 2(Xm _ Xr)/(3 /2-3x r )_ (56) 



1 + Z * \ C u. 

This expression was derived assuming 7 ^> 7d- Using the numerical values of Eq. fl2]), 
this corresponds to 

1 + " d(7) ^0.16 7 °- 2 . (57) 

The damping terms hence are more efficient for the larger length scales and the scaling 
solution takes over after an expansion factor of typically O(10), i.e. when Zd — 300. For 
7 C < 7 7d, one finds the same expression as in Eq. (1571) up to the replacement 7 — > y&. 
As a result, if 7d is very small, the scaling solution for the small loops may not overcome 
AC before today. For instance, if 7 < 7d = 10~ 14 , one finds z& ~ 0. Physically, it means 
that some of the small loops today may actually come from the radiation era. Notice 
that, in addition to the damping effect, there is the horizon cutoff N*{$. > 2t*) = 0. 
This imposes that all terms involving A/ - * exactly vanish for z < 2^(7) with 

1 + ^(7) _ f7 + 7d\ 2/3 
1 + "V2 + 7d/ 

However, compared to Eq. ( |57|) . z h < z d such that the horizon cutoff is less efficient 
than damping. Let us remark that for the typical value of GU — 10~ 7 , Eq. ([9]) implies 
that 7d — 10~ 5 and the residual loops from the radiation era disappear around z^ ~ 30. 
Finally, the second term involving a dependence on £*, the initial time for the matter 
epoch, in Eq. (|33|) evolves as 



(58) 



t*7 \a J V 1 + z, 



(59) 
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Figure 4. Relaxation of the loop number density distribution towards the matter era 
scaling from the radiation era scaling solution at z — . We have chosen the arbitrary 
values 7d = 10~ 4 , 7 C | rac i = 10~ 6 and 7 C i ma t = 10~ 7 while the other constants have 
been fixed to match numerical simulations [see Eq. J2J]. The relaxation takes longer 
for 7 ~ 7d and is completed at z ~ 70. Lower values of 7d would increase the relaxation 
time according to Eq. (|57"1) such that some small residual loops may still be present 
today. 

and contributes less than 10% as soon as z < z % ~ 65, independently of the value of 7d- 
In Fig. HI we have plotted the loop distribution in the matter era when TV* is given by 
the radiation era scaling solution. As expected from the above discussion, we find the 
relaxation time to be longer for loops with size 7 C < 7 ^ 7d while the redshift at which 
the matter scaling solution dominates is compatible with Eq. f )57|) . 

To conclude this section, the loop number density distribution today is expected 
to be in full scaling on all length scales provided 7d > 10 -14 . If 7<j < 10 -14 , then some 
of the loops having a size 7c < 7 ^ 7d have been formed in the radiation era and 
the distribution is accordingly modified. On the other length scales, the distribution 
matches with the matter era attractor. For those cases, one can use the full expression 
given in Eq. ( |33i) to get the precise loop number density distribution, as plotted in Fig. HJ 

3.3. Cosmological applications 

3.3.1. Density parameter of cosmic string loops We may now use our knowledge of 
the loop number density distribution in the matter era to compute various quantities 
describing the cosmological impact of these loops today For simplicity, we neglect the 
recent change in the universal expansion from matter to vacuum energy domination and 
continue to use Eq. (JSJ) to relate cosmic time to the redshift. Furthermore, we assume 
that 7d > 10~ 14 such that the distribution is in its (matter era) full scaling regime. From 



Cosmic string loop distribution on all length scales and at any redshift 



19 



Eqs. 032]) and ( 15"5|) . the loop number density distribution then simplifies to 

+4x7 >* ^ C A 1 7d 

t Fin > 7«*) 



7 2_2x (7 + 7d)V 2-2 X 7 + 7d 

^(7<7ct) = g^( 1+ 1 * W^^V, (60) 
7 + 7d V 2xc-27 + 7d/ \7 + 7d/ 

where C c and X are explicitly given in Eqs. f )25|) and f )29|) by using Eq. (155]) . and we 

no longer distinguish between x R and x M as only the latter is relevant for the following 

considerations. The energy density in the form of cosmic string loops of all sizes is given 

by 

Po = 5 I t 4 ^(7,*)7d 7 , (61) 
t Jo 

where the upper value comes from the distance to the horizon in the matter era using 
v — 2/3 (note that this expression does not include the Kibble's loops whose contribution 
remains, however, negligible). Although Eq. ( )6T|) can be explicitly integrated in terms 
of hypergeometric functions, the resulting expression is not particularly illuminating. A 
very good approximation can be obtained by piecewise integrating the asymptotic power 
law expansions of Eqs. fl36|) . f )38|) and f l4T|) in the three domains 7 < 7 C , 7c < 7 < 7d, 
and 7 > 7d- After some algebra, one finds 



Pot 2 C 



l-2x flc 



_A X il -2 X ) 4 X V7d 

7d xi-2x 



2x 



(62) 



(63) 



1 - 2 X V 3 

which can be further approximated for 7 C 7d <C 1 by 

Pot 2 C 1_ 

U ~ Ax{l-x){l-^x)l\~ 2x ' 
This expression clearly shows that the energy density in the form of loops is essentially 
fixed by the gravitational decay scale 7d and the scaling function properties through C 
and x- Smaller 7d implies larger values for p /U. Notice also that the scaling evolution 
is similar to that of the long string ones; the quantity t 2 p Q /U is time independent 
and universal. Doing the integration exactly and then expanding the hypergeometric 
functions gives the same result, up to the first numerical coefficient which changes by 
fifteen percents and reads CV/[(2 — 2%) sin(27rx)]. In terms of the density parameter 
Q = Po/Pcrit, where p cri t is the critical density for a flat universe today, we find 

(l- X )8in(27r X ) 7 J- 2 *' (64) 
which is valid at all scaling times during the matter era. Up to some small additional 
dilution factor coming from the recent cosmological constant domination, this expression 
may be extrapolated until today. Using Eq. iQ with F ~ 10 2 and the numbers derived 
before, one finds 

Q ^ 0.10 x [GU)°- 59 . (65) 
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This cannot exceed 10~ 5 with the values of GU currently allowed by CMB observations, 
i.e. GU < 7 x 10~ 7 , and hence Nambu-Goto cosmic string loops cannot provide a viable 
Dark Matter candidate, a scenario which has been dicussed in Ref. [57] . 

3.3.2. Number density of loops in a box of size L Again assuming the full scaling 
solution in the matter era, we can estimate the number density of cosmic string loops 
present around us today. Taking a box of size L, Eq. (1601) can be used to determine the 
number density of loops having a length i < L at t: 

rL/t 

t 3 n L = / t 4 J"(7,t)d 7 . (66) 
Jo 

Assuming the box is big enough, with 7c < 7d <C L/t, we find the dominant part of the 
loop density to be independent of L, 
3 2 Xc - 2 X C 

t n L = -7, 1 ■ 67 

2%c - 1 7d7 c 

As expected, the dependence in Xc is very weak such that the first term in the previous 
equation can be approximated to unity. Using Eqs. and (TT4"j) with T ~ 10 2 as before 
and T ~ 10, yields 

n L * 6 -^l {GU) -^. (68) 

Using t = to from Eq. (jSJ) evaluated for z = (which again amounts to neglecting the 
recent change to a vacuum energy dominated universe) with GU ~ 7 x 10~ 7 gives the 
lower bound 

n L ~ 5.5 x lO^Mpc" 3 . (69) 

This number is quite small, but interestingly it scales as a negative power of GU, such 
that loops from very light cosmic strings are the more numerous. Our result suggests that 
lensing events [68] from loops are more frequent for low values of GU (but also of smaller 
amplitude). Note that with the above string tension, the gravitational backreaction 
length is at 7c to — 8 pc, and the gravitational decay length at 7 dto — 380 kpc. 



4. Conclusions 



In this paper, we have proposed a phenomenological model of Nambu-Goto cosmic 
string loop evolution based on the Polchinski-Rocha loop production function, which 
we adjusted to fit numerical simulations. We have considered both loop decay by 
gravitational wave emission and string smoothing due to gravitational backreaction. We 
have found that the loop number density distribution rapidly assumes a universal form 
on all length scales, which ends up being insensitive to the details of the backreaction 
effects. It solely depends on the scale 7c below which these effects dominate. In the 
matter era, the energy density parameter of loops Q is a positive power of the string 
tension GU and ends up being completely negligible. Hence, cosmic string loops are 
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unfit to account for the Dark Matter contribution to our Universe's energy density. 
On the other hand, the density number of loops is boosted by a negative power of GU, 
meaning that very low tension loops could actually be abundant in the present Universe. 

Our work may be extended along several lines. For instance, since the universal 
loop distribution plotted in Fig. [2] is significantly different than the loop distributions 
used so far in the literature, various constraints on GU coming from gravitational wave 
emission, lensing or pulsar-timing should be re-examined [321 [331 Ell |3H [70]. All along 
the calculation, we have kept 7d and 7 C as undetermined parameters (apart from the 
numerical examples) since these parameters may equally well be used to represent other 
physical effects. For instance, 7a does not need to be given by Eq. (jHJ) but could represent 
the decay rate of another energy loss mechanism, such as particle emission, provided it 
can still be described phenomenologically by Eq. (JSJ). Another extension would be the 
implementation of a low intercommutation probability (as in the PR model, we have set 
this probability P = 1 throughout this work) to extend our approach to the evolution 
of cosmic superstring loops. 



Acknowledgments 

It is a pleasure to thank Patrick Peter and Daniele Steer for enlightening discussions. 
The work of L. L. and C. R. is partially supported by the Belgian Federal Office for 
Science, Technical and Cultural Affairs, under the Inter-university Attraction Pole 
Grant No. P6/11. The work of M. S. is partially supported by the Sciences & 
Technology Facilities Council (STFC-UK), Particle Physics Division, under the grant 
ST/G000476/1 "Branes, Strings and Defects in Cosmology". 



References 

T. W. B. Kibble, Topology of cosmic domains and strings., J. Phys. A 9 (1976) 1387-1398. 
M. B. Hindmarsh and T. W. B. Kibble, Cosmic strings, Rept. Prog. Phys. 58 (1995) 477-562, 



hep-ph/9411342 



A. Vilenkin and E. P. S. Shellard, Cosmic Strings and Other Topological Defects. Cambridge 

University Press, Cambridge, England, July, 2000. 
M. Sakellariadou, Cosmic strings, Led. Notes Phys. 718 (2007) 247-288, |hep-th/0602276| . 



D. Kirzhnits and A. Linde, Macroscopic consequences of the Weinberg model, Phys. Lett. B 42 
(Dec., 1972) 471-474. 

R. Jeannerot, J. Rocher, and M. Sakellariadou, How generic is cosmic string formation in SUSY 

GUTs, Phys. Rev. D68 (2003) 103514, |hep-ph/0308134| . 
J. Rocher and M. Sakellariadou, D-term inflation, cosmic strings, and consistency with cosmic 

microwave background measurement, Phys. Rev. Lett. 94 (2005) 011303, hep-ph/041 2143| . 
C. P. Burgess et. at, The Inflationary Brane-Antibrane Universe, JEEP 07 (2001) 047, 

|hep-th/0105204| . 

S. Sarangi and S. H. H. Tye, Cosmic string production towards the end of brane inflation, Phys. 

Lett. B536 (2002) 185-192, |hep-th/0204074| . 
G. Dvali and A. Vilenkin, Formation and evolution of cosmic D-strings, JCAP 0403 (2004) 010, 



hep-th/0312007 . 



Cosmic string loop distribution on all length scales and at any redshift 



22 



N. T. Jones, H. Stoica, and S. H. H. Tye, The production, spectrum and evolution of cosmic 
strings in brane inflation, Phys. Lett. B563 (2003) 6-14, |hep-th/0303269|. 

A.-C. Davis and T. Kibble, Fundamental cosmic strings, Contemp. Phys. 46 (Sept., 2005) 
313-322, |hep-th/0505050| . 

M. Sakellariadou, Cosmic Superstrings, Phil. Trans. Roy. Soc. Lond. A366 (2008) 2881-2894, 
[arXiv: 0802. 3379] . 

E. J. Copeland and T. W. B. Kibble, Cosmic Strings and Superstrings, Proc. Roy. Soc. Lond. 

A466 (2010) 623-657, [arXiv:0911 . 1345] . 
M. Sakellariadou, Cosmic Strings and Cosmic Superstrings, Nucl. Phys. Proc. Suppl. 192-193 

(2009) 68-90, [arXiv: 0902 .0569] . 

A. Albrecht and N. Turok, Evolution of cosmic string networks, Phys. Rev. D40 (Aug., 1989) 
973-1001. 

D. P. Bennett and F. R. Bouchet, Cosmic-string evolution, Phys. Rev. Lett. 63 (Dec., 1989) 
2776-2779. 

D. P. Bennett and F. R. Bouchet, High- resolution simulations of cosmic-string evolution. I. 
Network evolution, Phys. Rev. D41 (Apr., 1990) 2408-2433. 

B. Allen and P. Shellard, Cosmic-string evolution - A numerical simulation, Phys. Rev. Lett. 64 
(Jan., 1990) 119-122. 

N. Jarosik et. at, Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 

Sky Maps, Systematic Errors, and Basic Results, arXiv : 1001 . 4744, 
N. Bevis, M. Hindmarsh, M. Kunz, and J. Urrcstilla, Fitting CMB data with cosmic strings and 

inflation, Phys. Rev. Lett. 100 (2008) 021301, |astro-ph/0702223| . 

A. A. Fraisse, C. Ringeval, D. N. Spergel, and F. R. Bouchet, Small-Angle CMB Temperature 
Anisotropics Lnduced by Cosmic Strings, Phys. Rev. D78 (2008) 043535, [arXiv : 0708 . 1162] . 

K. Takahashi et. ai, Non-Gaussianity in Cosmic Microwave Background Temperature 

Fluctuations from Cosmic (Super-) Strings, JCAP 0910 (2009) 003, |arXiv:08irT 4698j. 
M. Hindmarsh, C. Ringeval, and T. Suyama, The CMB temperature bispectrum induced by 

cosmic strings, Phys. Rev. D80 (2009) 083501, [arXiv: 0908. 0432] . 
M. Hindmarsh, C. Ringeval, and T. Suyama, The CMB temperature trispectrum of cosmic 

strings, Phys. Rev. D81 (2010) 063505, [arXiv : 0911 . 1241] . 
D. M. Regan and E. P. S. Shellard, Cosmic String Power Spectrum, Bispectrum and Trispectrum, 

Phys. Rev. D82 (2010) 063527, [arXiv: 0911 .2491] . 
D. Yamauchi et. ai, Skewness in CMB temperature fluctuations from curved cosmic 

(super-) strings, arXiv: 1004.0600. 
M. Landriau and E. P. S. Shellard, Cosmic String Lnduced CMB Maps. llarXiv : 1004 . 2885. 
N. Bevis, M. Hindmarsh, M. Kunz, and J. Urrestilla, CMB power spectra from cosmic strings: 

predictions for the Planck satellite and beyond, Phys. Rev. D82 (2010) 065004, 

[arXiv: 1005.2663] . 

C. Ringeval, Cosmic strings and their induced non-Gaussianities in the cosmic microwave 
background, Adv. Astron. 2010 (2010) 380507, [arXiv: 1005 .4842] . 

B. Allen and E. P. S. Shellard, Gravitational radiation from cosmic strings, Phys. Rev. D 45 
(Mar, 1992) 1898-1912. 

T. Damour and A. Vilenkin, Gravitational wave bursts from cusps and kinks on cosmic strings, 

Phys. Rev. D64 (2001) 064008, |gr-qc/0 104026] . 
X. Siemens et. ai, Gravitational wave bursts from cosmic (super) strings: Quantitative analysis 



and constraints, Phys. Rev. D73 (2006) 105001, |gr-qc/0603115| . 
S. Olmez, V. Mandic, and X. Siemens, Gravitational-Wave Stochastic Background from Kinks 
and Cusps on Cosmic Strings, arXiv : 1004 . 0890 

F. Dubath and J. V. Rocha, Periodic gravitational waves from small cosmic string loops, Phys. 
Rev. D76 (2007) 024001, |gr-qc/0703"l09] . 

G. Vincent, N. D. Antunes, and M. Hindmarsh, Numerical Simulations of String Networks in the 



Cosmic string loop distribution on all length scales and at any redshift 



23 



Abelian-Higgs Model, Phys. Rev. Lett. 80 (Mar., 1998) 2277-2280, [hep-ph/9708427| . 
[37] J. N. Moore, E. P. S. Shellard, and C. J. A. P. Martins, Evolution of Abelian-Higgs string 

networks, Phys. Rev. D65 (Jan., 2001) 023503, |hep-ph/0107171| . 
[38] M. Hindmarsh, S. Stuckey, and N. Bevis, Abelian Higgs Cosmic Strings: Small Scale Structure 

and Loops, Phys. Rev. D79 (2009) 123504, |arXiv:0812. 1929j . 
[39] G. R. Vincent, M. Hindmarsh, and M. Sakellariadou, Scaling and small scale structure in cosmic 

string networks, Phys. Rev. D56 (1997) 637-646, |astro-ph/9612135| . 
[40] V. Vanchurin, K. Olum, and A. Vilenkin, Cosmic string scaling in flat space, Phys. Rev. D72 

(2005) 063514, |gr-qc/0501040| . 

[41] C. Ringeval, M. Sakellariadou, and F. Bouchet, Cosmological evolution of cosmic string loops, 

JCAP 0702 (2007) 023, | astro-ph/0511646| . 
[42] V. Vanchurin, K. D. Olum, and A. Vilenkin, Scaling of cosmic string loops, Phys. Rev. D74 

(2006) 063527, |gr-qc/0511159] . 

[43] C. J. A. P. Martins and E. P. S. Shellard, Fractal properties and small-scale structure of cosmic 

string networks, Phys. Rev. D73 (2006) 043515, |astro-ph/0511792| . 
[44] K. D. Olum and V. Vanchurin, Cosmic string loops in the expanding universe, Phys. Rev. D75 

(2007) 063521, astro-ph/0610419| . 

[45] A. Vilenkin, Gravitation radiation from cosmic strings., Phys. Lett. B 107 (1981) 47-50. 
[46] M. Sakellariadou, Gravitational waves emitted from infinite strings, Phys. Rev. D42 (1990) 
354-360. 

[47] M. Hindmarsh, Gravitational radiation from kinky infinite strings, Phys. Lett. B251 (1990) 
28-33. 

[48] X. Siemens and K. D. Olum, Gravitational radiation and the small-scale structure of cosmic 

strings, Nucl. Phys. B611 (2001) 125-145, | gr-qc/0104085| . 
[49] X. Siemens, K. D. Olum, and A. Vilenkin, On the size of the smallest scales in cosmic string 

networks, Phys. Rev. D66 (2002) 043501, |gr-qc/0203006l . 
[50] J. Polchinski and J. V. Rocha, Cosmic string structure at the gravitational radiation scale, Phys. 

Rev. D75 (2007) 123503, gr-qc/0702"055] . 
[51] D. Austin, E. J. Copeland, and T. W. B. Kibble, Evolution of cosmic string configurations, Phys. 

Rev. D48 (1993) 5594-5627, |hep-ph/9307325| . 
[52] C. J. A. P. Martins and E. P. S. Shellard, Quantitative String Evolution, Phys. Rev. D54 (1996) 

2535-2556, |hep-ph/9602271| . 
[53] E. J. Copeland, T. W. B. Kibble, and D. A. Steer, The evolution of a network of cosmic string 

loops, Phys. Rev. D58 (1998) 043508, |hep-ph/9803414| . 
[54] C. J. A. P. Martins and E. P. S. Shellard, Extending the velocity-dependent one-scale string 

evolution model, Phys. Rev. D65 (2002) 043514, |hep-ph/0003298| . 
[55] J. Polchinski and J. V. Rocha, Analytic study of small scale structure on cosmic strings, Phys. 

Rev. D74 (2006) 083504, |hep-ph/0606205| . 
[56] F. Dubath, J. Polchinski, and J. V. Rocha, Cosmic String Loops, Large and Small, Phys. Rev. 

D77 (2008) 123528, [arXiv: 071 1.0994] . 
[57] V. Vanchurin, Non-linear dynamics of cosmic strings with non-scaling loops, Phys. Rev. D82 

(2010) 063503, [arXiv: 1005. 1680] . 
[58] V. Vanchurin, Semi-scaling cosmic strings, arXiv : 1005 . 5721. 

[59] J. V. Rocha, Scaling solution for small cosmic string loops, Phys. Rev. Lett. 100 (2008) 071601, 
[arXiv: 0709.3284] . 

[60] E. Komatsu et. al, Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 

Cosmological Interpretation, arXiv: 1001 .4538. 
[61] E. J. Copeland and T. W. B. Kibble, Kinks and small-scale structure on cosmic strings, Phys. 

Rev. D80 (2009) 123523, [arXiv : 0909. I960] . 
[62] T. Vachaspati and A. Vilenkin, Formation and Evolution of Cosmic Strings, Phys. Rev. D30 

(1984) 2036. 



Cosmic string loop distribution on all length scales and at any redshift 



24 



[63] G. Karra and R. J. Rivers, Initial vortex densities after a temperature quench, Phys. Lett. B414 

(1997) 28-33, |hep-ph/9705243| . 
[64] M. Hindmarsh and A. Rajantie, Phase transition dynamics in the hot Abelian Higgs model, Phys. 

Rev. D64 (2001) 065016, |hep-ph/0103311| . 
[65] T. Vachaspati and A. Vilenkin, Formation and evolution of cosmic strings, Phys. Rev. D30 

(Nov., 1984) 2036-2045. 

[66] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products. Academic Press, 

New York and London, 1965. 
[67] Y. Cui and D. E. Morrissey, Non-Thermal Dark Matter from Cosmic Strings, Phys. Rev. D79 

(2009) 083532, [arXiv: 0805. 1060] . 
[68] A. A. de Laix and T. Vachaspati, Gravitational lensing by cosmic string loops, Phys. Rev. D54 

(1996) 4780-4791, |astro-ph/9605171| . 
[69] M. S. Pshirkov and A. V. Tuntsov, Local constraints on cosmic string loops from photometry and 

pulsar timing, Phys. Rev. D81 (2010) 083519, |arXiv: 0911 . 4955j . 
[70] A. V. Tuntsov and M. S. Pshirkov, Quasar variability limits on cosmological density of cosmic 

strings, Phys. Rev. D81 (2010) 063523, [arXiv: 1001 .4580] . 



